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Abstract 

A tomographic technique has been used in the past to decompose 
complex signals in its components. The technique is based on spectral 
decomposition and projection on the eigenvectors of a family of unitary 
operators. Here this technique is also shown to be appropriate to 
obtain the instantaneous phase derivative of the signal components. 
The method is illustrated on simulated data and on data obtained 
from plasma reflectometry experiments in the Tore Supra. 
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1 Introduction: Plasma density from reflec- 
tometry and its multicomponent nature 

Density measurements play an important role in the study and operation 
of magnetically confined plasmas. Microwave reflectometry is a radar-like 
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technique which infers the plasma density from the reflection on the (cutoff) 
layers where the refractive index vanishes pQ. For example, for propagation 
perpendicular to the magnetic field with the electric field of the wave parallel 
to the magnetic field in the plasma (0-mode), the refractive index is 



^' V Bom, {27Tff ^ ' 

where is the electron density, uOp = ^^^^^7^ j ^ the plasma frequency, e 
and rrie the electronic charge and mass, Eq the permittivity of the vacuum 
and / = ^ the frequency of the probing wave. When the plasma frequency 
equals the probing frequency the index of refraction vanishes, the wave is 
reflected and the density ric of the cutoff layer may be derived from 

Uc = (2) 

Mixing the reflected wave Ef{{t) with the (reference) incident wave EQ{t), 
the mixer output is 

^ {E^, {t) + El (t)) + Eo (t) Er (t) cos (t) 

In the interference term Eq (t) Er (t) cos (t) , Eq (t) Er (t) depends on many 
factors, microwave generator power, plasma scattering properties, turbu- 
lence, etc., therefore it is (t) that contains the most reliable information 
about the plasma density. 

The location Xc (ujp) of the reflecting layer for the frequency Up is related 
to the group delay 

_ d(j) (uj) _ l_d^ 

dw ~ 271 df ^ ^ 



by (0-mode) 



Xc [ujp) = Xo + - / duj -— (4) 



For a linear frequency sweep of the incident wave 

fit) = fo + lt (5) 
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one obtains 



d(j) 1 d(f) 
df 7 dt 



f 



(6) 



Therefore, measurement of the plasma density hinges on an accurate deter- 
mination of the "instantaneous frequency" Accuracy in the measurement 
of this quantity is quite critical because, the location of the reflecting layer 
being obtained from the integral in (jl]), errors tend to accumulate. 

Several methods have been devised to obtain the group delay Tg from the 
refiectometry data (for a review see ^j). Among them, time- frequency anal- 
ysis [3j has been, so far, the most promising technique. The Wigner-Ville 
(WV) distribution [1] [5] although providing a complete description of the 
signal in the time-frequency plane, raises difficult interpretation problems 
due to the presence of many interference terms that impair the readabil- 
ity of the distribution. This occurs because the VW (quasi-) distribution is 
not a probability distribution, has complex amplitudes and may have large 
amplitude values in frequency regions which are not contained in the signal 
spectrum. For this reason the time-frequency method that has been pre- 
ferred is the spectrogram [3] [6] [7j [8], that is, the squared modulus of the 
short-time Fourier transform 



h (u) being a peaked short-time window. 

The spectrogram does not really provide the instantaneous frequency, 
because that notion is not well defined anyway. All it gives is the product 
of the spectra of x (t) and h {t) . The way the spectrogram is used to infer 
the local rate of phase variation ^ is to identify this quantity with the 
maximum or the with the first moment of the spectrogram. An additional 
problem comes about because unwanted phase contributions due to plasma 
turbulence may have a higher amplitude than the contributions due to the 
profile. Correction techniques have been developed to compensate for this 
errors, based for example on Floyd's best path algorithm. The choice of the 
window function is also an important issue and, in particular, an adaptive 
spectrogram technique has been developed to maximize the time-frequency 
concentration [2]. 

In addition to the delicate nature of the extraction of the phase derivative 
from an interference signal, another important question is the multicompo- 
nent structure of this signal. The signal that is actually received contains, in 
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SP (t, /) = / x{u)h{u~t)e 
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addition to the reflection on the plasma, reflections on the porthole and multi 
reflections of the waves on the wall of the vacuum vessel. Separation of these 
latter components from the plasma reflections is an essential step to obtain 
reliable density results. Separation by frequency flltering is not appropri- 
ate because there is considerable frequency overlap between the components. 
In a previous paper [H] we have developed a method to separate the signal 
components based on a tomographic representation[Tn] [H], which gives a 
positive density Mf{x,9) of the signal along all possible 6'-directions in the 
time-frequency plane. 

The tomogram representation Mf{x,6) gives, for ^ = the time rep- 
resentation of the signal, fit), and for 9 = ^ the frequency representation, 

/(z/). Frequency flltering corresponds to component separation of the sig- 
nal at 6' = I and, from the many examples that were studied, one concludes 
that, in general, this is not the most convenient direction to isolate the signal 
components. For example, for the reflectometry signals that were studied, 
we have more information if the separation of the components is performed 
at 6^ = than at 6* = |. 

In the next section we flrst make a brief review of the tomographic method 
for component separation and then, using the same mathematical framework, 
show how one can obtain the phase derivative from the isolated components. 

Finally, in the last two sections, the methods are applied both to simu- 
lated data and to actual reflectometry data collected in the Tore Supra. 



2 Tomograms, components and the phase deriva- 
tive 

In [9] we described in much detail the use of tomograms for the component 
factorization of complex signals. Here we just recall some basic facts for the 
reader's convenience. 

We deflne a (time-frequencj|l]) tomogram as a family of probability dis- 
tributions, Mf{x,9) associated to any signal fit), t e [0, T] by 



Msix,e) 



fim7it)dt 



l</,*:'^>l 



(7) 



^As explained in other non-commuting operator pairs may be chosen 
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with 

^'fit) = ^exp{—-—t^ + —-t] 8 
^/T V 2 sm 6* sm9 J ^ ^ 

Notice that >< ^^'"^| are spectral projections of an unitary operator 

U (6) and therefore ([7]) performs a spectral decomposition of the signal. 

First we select a subset of numbers in such a way that the corre- 
sponding family |\E'^'J(t)}^ is orthogonal and normalized: 

<<'L^l':>=Sm,n (9) 

A glance at the shape of the functions ([8]) shows that, for fixed 6, the 
oscillation length at a given t decreases when |x| increase. As a result, the 
projection of the signal on the {\Ef^^(t) } basis locally explores different scales. 
On the other hand the local time scale is larger when 6 also becomes larger, 
in agreement with the uncertainty principle for a non-commuting pair of 
operators. 

We then consider the projections of a signal f{t) 

4(/)=</'^'f > (10) 

and use the coefficients c^xAf) signal processing purposes. 

In particular, a multi-component analysis of the signal |2| is done by select- 
ing subsets J-'k of the and reconstructing (/c-component) partial signals 
by restricting the sum to 

m = E (11) 

for each fc-component. 

In the present work we analyze the phase derivative of a complex signal 
f{t) = A{t)e'^'^^^^ and consider the cases where f{t) already corresponds to 
one of the components determined as in [9j. That is, after a convenient 
factorization of the signal is performed, the search for the phase derivative is 
made on each component. 

In the reflectometry technique the experimental signal is already complex 
(it consists of one recorded interference term composed of in-phase and 90° 
phase shifted signals). Therefore we have no ambiguity in the definition of 
the amplitude and phase of f{t). For other type of signals, where only the 
real part is available, the construction of a complementary imaginary part 
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is an usual technique for which there are standard methods available in the 
signal analysis literature (see [T2| and references therein). 

Given a signal f{t) = A(t)e''^^^^ the time derivative of the phase may be 
obtained from 



d 



df 



^^W = /m|j^l (12) 



For our decomposed components one has 



^ <P{t) = Im(l^] (13) 



with 



T(t) = $:4(/)|v,^f(t) (14) 



and 

^w = E4(/)*^.fw (15) 

Notice that an explicit analytic expression for ^\l/^'J(t) is known, namely: 

^^l'^{t) = I ( + ^'f{t) (16) 

(9t ^" ^ ^ V sin^ sin^y ^" ^ ^ ^ ' 

and therefore we obtain a direct expression for the phase derivative in terms 
of the coefficients c^^S^) without having to use the values of / for neighboring 
values of t. This provides a more robust method to estimate the derivative. 
We call the Tomographic Direct Method (TDM) the method of the 
computation of the phase derivative of f{t) using f[T^ . 

Notice that in the calculation of the imaginary part in f[T^ the value of 
the amplitude A (t) plays no role. Therefore we may use what will be called 
a Tomographic Normalized Method (TNM) defined in the same way 
as TDM but with a normalized signal replacing f{t). For calculations 
on the signal carried out with absolute precision the results of TDM and 
TNM should coincide. However because of numerical errors, normalization 
of the signal amplitude, before further processing, might have some merit 
mostly in the small amplitude regions. 

There are still two specific issues to be addressed when dealing with the 
reconstruction of the phase derivative of f{t). The first is a general problem 
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in signal analysis, namely denoising. We recall from [H] that Tomogram- 
Based Denoising (TBD) consists in eliminating from (1151) the c^.^(/) such 
that 

i4(/)r<^ (17) 

for some chosen threshold e. 

Another way, often used for denoising, consists in locally smoothing the 
signal by computing a Local Mean (LMm) G of a function g{t), known in 
the signal processing community as moving average FIR filter of order 
2m + 1 by: 

k=—m 

The second issue is how to handle the difficult problem of accurate mea- 
surement of the phase, hence also of the phase derivative, when the signal 
amplitude is very small. Given a complex signal f{t) we define the trun- 
cated Phase Derivative (tPD) by 

if /(t) < a, for some convenient threshold a, 
else 

^%{t) = Im(^] (19) 



5t Km 

Notice that tPD simply sets the value of the phase equal to a constant when 
the signal amplitude prevents its accurate estimation. 

In the following sections we present some advantages and drawbacks of 
these tools by applying them to several simulated and experimental signals. 



3 Examples: Simulated data 

In this section we apply the general method to two types of simulated sig- 
nals. The first example shows how the phase derivative of a sinusoidal signal 
may be computed with accuracy, even when noise is present. In the second 
example, we focus on the phase derivative of signals with non linear phase. 
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Our data consists of complex functions y{t) = A(t)e*'^^*'' with phase and 
phase derivative dt4>{t) unambiguously defined. The analysis of all the sim- 
ulated signals is based on tomograms with ^ = | as for the same data in 

For the simplest case, the signal is: 

y{t) = exp{i75t),t e [0, 20] s (20) 

TDM alone gives an excellent result (mean value of the computed dt4>{t) is 
74:.Srd/s and the standard deviation (sdev) 0.3rd/ s to be compared with the 
Fourier Transform for which the resolution is equal to A/ = ~ O.Slrd/s. 

If we add a (complex) noise b{t) to with SNR = 10d£§, TDM, 
not surprisingly, still shows a good mean result (75.9 rd/s) but has a large 
uncertainty (sdev=40 rd/s). The use of LM alone is not sufficient in this case 
(sdev=3.5 rd/s for a LM15) but TBD allows a TDM with great accuracy 
(sdev=0.8 rd/s) that may even be improved by a ultimate use of a LM, 
(sdev=0.6 rd/s for LM5). 

It is worthwhile to mention how denoising using the spectral decomposi- 
tion of the operator U{a) (TBD) works so efficiently, a result that is also 
confirmed in the subsequent examples. 

We proceed to the analysis of a signal which aims to mimic, in a simplified 
way, the case of an incident plus a refiected wave delayed in time and with 
an acquired time-dependent change in the phase. In this case the simulated 
signal y{t) is the sum of an "incident" chirp yo{t) and a "deformed refiected" 
chirp ynit). Noise is added to the signal and the SNR = lOdB. However 
thanks to the analysis in [9] we may consider these two waves separately. 

For the "incident" chirp yo{t) the analysis is performed in two different 
situations that differ mainly in an amplitude term. 

The signal is: 



yo{t) = b{t),te[0,3]s 

y,{t) = A(t)e^*»W+6(t),te [3,18]s 

yoit) = 6(t),te [18,20] s (21) 



2The SNR is defined by : SNR{y,b) = lOlogio^ with Py ^ \y{t)\^ dt and 
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where ^o(t) = (^ot"^ + bot and oq, &o are chosen to have dt4>{3) = 75rd/s and 
at0(18)=5O rd/s. 

Here A{t) is one in the first case and in the second case, A{t) defined by 
fl22|) is defined for t E [3, 18] s by equation fl22l) and presented on Fig{Tl Here 

= 6 and Uk is randomly chosen between rd/s and 7.5 rd/s. Notice that 
for t=14 s, A is very small. 




0.3 - 
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Figure 1: Amplitude A(t), defined by equation (j22|) . of the signal y{t) defined 
by equation fl2T]) . 

For this signal, tPD considerably improves the result for t G [0, 3] and 
t G [18, 20], as it is easy to understand since the phase derivative of a random 
signal may have large local values but the corresponding amplitude of the 
total signal be small. 

After using the tPD, we summarize the performances of the different 
tools in the following Table 1 in terms of their S(ietl§. 

In Figl2] we show the graphic representation of the reconstructed phase 
derivative, for A{t) ^ 1, corresponding to TDM+LM (LM5) and TDM+TBD 

■'in this case the standard deviation sdev is defined using the difference between the 
analytic expression of the known dt(j){t) and the corresponding estimated value 
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sdev 


TDM 


TL)M + LMs 


TL)M + 


TDM + TfiD + LM5 


A{t) = 1 


38.5 rd/s 


4.5 rd/s 


1.8 rd/s 


1.5 rd/s 


Ait) ^ 1 


51 rd/s 


11.2 rd/s 


1.9 rd/s 


1.5 rd/s 




TNM 


TNM + LM5 


TA^M + TBD 


TNM + TBD + LM5 


A{t) = 1 


23.5 rd/s 


3.9 rd/s 


1.8 rd/s 


1.5 rd/s 


Ait) ^ 1 


39.3 rd/s 


10.7 rd/s 


1.3 rd/s 


2.2 rd/s 



Table 1: Comparison of the different tools in terms of their sdev for the 
signal yo defined by fl2T|) . 

of Table 1. As can be seen, the combined use of the tools described in section 
2 allows a very efficient reconstruction of the phase derivative in this case, 
except when the signal is very small, for t ~ 14s. In particular TDM (or 
TNM )+ TBD gives very good results for an amplitude varying signal. It is 
however worthwhile to mention that the tomogram spectral family ([HD is par- 
ticularly well adapted to this type of "incident wave" since in the limit case 
of an infinite time domain the corresponding spectrum would be reduced to 
a unique c^(/). But if, on one hand, we take advantage of this fact because 
the incident wave in reflectometry has this shape, on the other hand, the 
next example shows that the good performances of the tool are not limited 
to this particular non linear phase shape. 

We also notice from Table 1 that, even before filtering, the normalisation 
TNM improves the results. This arises mostly from the processing of the 
small amplitudes regions. 

Let now consider the "deformed reflected chirp" ynit) defined by : 

ynit) = A(t)e**«« + t G [0, 20] s (23) 

where $_R(t) = aut'^ + h^t + lOti and a^j, hji are chosen to have (9j0(O) = 
75rd/s and 9(0(2O)=5O rd/s. As before A{t) is one in the first case and 
defined by (1221) in the second case with, in each case, a noisy component b{t) 
with SNR = 1MB. 

The performances of the different tools, summarised in Table 2 in terms 
of their sdev, show how the tomogram based tools perform very accurate 
estimations of the (local) phase derivative in cases where other methods may 
have some difficulties. In particular, TBD seems a very efficient method to 
denoise the signal as it can be seen in the Figj3l 
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Figure 2: Phase derivative of yo{t), defined by equation fl2T]) for the case 
A{t) ^ 1, using TDM+LM (LM5) and TDM+TBD (bold hue) on the 
tomogram for ^ = |. 

4 Application to reflectometry data 

We now show the abihty of the tomographic methods to extract the phase 
derivative of an experimental signal coming from reflectometry measurements 
during a discharge in the Tore Supra at Cadarache. 

The sweep- frequency reflectometry system of Tore Supra launches a prob- 
ing wave on the extraordinary mode polarization (X mode) in the V band 
(50-75 GHz) [6], [7], [13]. The emitting and receiving antennas are located 
at about 1.20 m from the plasma edge, outside the vacuum vessel. The reflec- 
tometry system repeatedly sends sweeps of duration 20/is. The heterodyne 
reflectometers, with I jQ detection, provide a good Signal to Noise Ratio, up 
to 40(i5. For each sweep, the reflected chirp E^it) is mixed with the incident 
sweep -Eo(^) and only the interference term is recorded as an in-phase and a 
90° phase shifted sampled signals: 

xi(t) = AQAn{t) cos(v9(t)) 
X2(t) = AQAn{t) sin((^(t)) 
Let the reflected signal be 

yit) = x,{t) + tX2{t) = A(t)e^'^W (24) 
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sdev 


TDM 


TL)M + LMs 


TBD 


+ LM5 


A{t) = 1 


27.1 rd/s 


6.6 rd/s 


2.0 rd/s 


1.5 rd/s 


Ait) ^ 1 


92 rd/s 


24.7 rd/s 


3.0 rd/s 


2.9 rd/s 




TNM 


TNM + LM5 


TA^M + TED 


TA^M + + LM5 


A{t) = 1 


27.1 rd/s 


4.7 rd/s 


2.0 rd/s 


1.5 rd/s 


Ait) ^ 1 


36.1 rd/s 


14.1 rd/s 


2.2 rd/s 


1.4 rd/s 



Table 2: Comparison of the different tools in terms of their sdev for the signal 
Ur defined by fl23|) . 

The phase derivative of the signal corresponding to the plasma component 
of y{t) is used to localize the cut-off density in the plasma. The amplitude 
of this signal A{t) = AoAfiit) corresponds to a low frequency. The real part 
of the signal y{t) is shown in Fig HI 

The tomogram at 6 = ^ was used to perform the factorization of the 
signal in [9]. Cutting the spectrum at e = 0.05 max(|c^^(|/)|) the signal is 
factorized in three components as shown on FiglH 

4.1 First component, the reflection on the porthole 

The first component, yi{t), of the reflectometry signal is a low frequency 
signal corresponding to the heterodyne product of the probe signal with the 
reflection on the porthole [13]. This complex signal is written : 

yi(t) = Ai(t)e^^^W (25) 

The phase derivative dt^pi{t) may be positive and proportional to the 
time Ti of the reflection of the probe signal on the porthole. If not, the 
reflectometry signal y{t) defined by ( IMl) is multiplied by e*"* for some a to 
calibrate the measurement. The real part and the modulus of yi{t) are shown 
in Figini 

The phase derivative of yi{t) was then computed using TDM and shown 
in FigJTl The mean value of the phase derivative is equal to —16.4 rd/s and 
its sdev to 2Ard/s (less than 15 %) in agreement with a rough estimation 
based on a spectrogram technique (see section 5). Since the phase derivative 
of yi{t) is negative, the reflectometry signal y{t) had to be calibrated to set 
dt^iit) proportional to ri (see conclusion). 
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Figure 3: Phase derivative of yR(t), defined by equation ([23]) for the case 
A{t) ^ 1, using TDM+LM (LM5) and TDM+TBD (bold hue) on the 
tomogram for ^ = |. 

We also shift the phase derivative of the other components by the same 
value. 

4.2 Second component, the plasma signal 

The second component has a Fourier spectra that fits the expected behavior 
corresponding to the reflection of the wave inside the plasma of the Tore 
Supra [13]. This component, 2/2 (^) is defined as: 

7/2(t) = A2(t)e^^^« (26) 

The modulus and real part are displayed together in the same plot (FigJHj). 

Even if the modulus is of low frequency in comparison to the real part of 
y2^\ the modulus is not constant. In particular the signal is very small in 
the first half. 

The amplitude of the signal yi^ for t G [0, ll]s, is very small in comparison 
to the amplitude for t G [ll,20]s. The power ratio of the signal ?/2 for 
t G [0, ll]s to the one for t G [11, 20]s equals 0.0350 (SNR ^ -15 dB). Then 
the first part of the signal can be considered as filtered noise and , with tPD, 

''The power Pg of a signal sit) is defined by : 
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Figure 4: Time representation of the reflectometry signal (real part). 

the phase derivative will not be computed. This fact may correspond to the 
difficulty of the incident wave to reach the plasma in the first (lower) band of 
"instantaneous frequencies" and therefore to a bad accuracy in the detection 
of the low densities present at the border of the plasma in tokamaks. 

The phase derivative of the last part of the signal, for t G [11, 20]s, cor- 
responding to TDM and TDM+LM are shown in FigjOl 

For comparison, the LMi^ filtered phase derivative of the TDM and of 
the TNM (bold), for t G [ll,20]s, are plotted on the same Figlini The 
results are nearly the same, except for small differences for t ~ 19.5s where 
the amplitude of the signal is very small. 

We conclude that combining TDM (or TNM) with LM gives an ac- 
curate estimation of the phase derivative of the plasma component. TNM 
appears to be performant when the amplitude of the signal small. This claim 
will be confirmed by a comparison with the usual spectrogram analysis in 
section 5. 
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Figure 5: Spectrum c^^(?/) of the reflectometry signal y{t) for do = % used 
in the factorization. 

4.3 Third component, the multireflection 

The last component of the reflectometry signal corresponds, [13] [9], to mul- 
tireflections of the wave on the wall of the vacuum vessel. This component, 
Vsit) is written as : 



The modulus A^it) and the real part of ysit) are presented together on 
the same figure (Fig lTTl) . As compared to the real part of 1/3 (t), the modulus 
Aslt) is a low frequency signal. We notice that for t > 13s the modulus is 
very small. 

The phase derivative of 1/3 (t) estimated using TDM+LM is plotted in 
Figini 

4.4 Components comparison 

The LM filtered phase derivative of the three components of the reflectome- 
try signal are plotted together on the same figure (Fig JTSjl . It is instructive to 
compare these phase derivatives. For the first component, the phase deriva- 
tive dt(pi{t) is almost constant. This is because the phase (pi{t) is mainly due 
to a simple reflection on a nearby object, the porthole. The reflection on the 
plasma is quite complex. The first part of the signal should be considered as 



y,{t) = A,{t)e" 



(27) 
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Figure 6: First component of the reflectometry signal corresponding to the 
reflection on the porthole (modulus and real part). 

filtered noise and shows that there is a problem in reconstructing the den- 
sity profile corresponding to this part of the sweep. The phase derivative of 
the third component of the signal, corresponding to multi reflections on the 
vessel, presents some similarities with the phase derivative of the first compo- 
nent, except for t > 14.5s. It will eventually be interesting to factorize again 
this component if some information related to the properties of the plasma 
close to the vessel walls can be extracted from it. The modulus of the three 
components are low frequency signals, compared to the signals themselves. 
Usually, the phase derivative obtained by TDM is accurate when filtered by 
LM. In this case TBD does not seem adequate for denoising purposes, be- 
cause it correlates with the component analysis and may eventually corrupt 
the factorization of the signal. 

5 Tomograms and spectrogram analysis 

In this section we obtain the "frequency" of the signal as a function of time, 
obtained by a moving window FFT spectrogram, and compare it with the 
phase derivative obtained by the tomographic techniques described before. 
The spectrogram is computed with a 64 points length window (the grain) and 
a 75% overlap rate using the maximum pick method [T2], allowing a FFT res- 
olution of lOrd/s on 121 time points. We avoided an estimation with higher 
time resolution because of FFT resolution constraints. For the tomographic 



16 



O 5 10 15 20 

time 



Figure 7: Phase derivative of the first component of the refiectometry signal 
corresponding to the reflection on the porthole, using TDM on the tomogram 
for 6* = — 

phase derivative estimation we used, as usually, TDM (or TNM) together 
with LM filtering. 

For the simulated "deformed reflected chirp" yR{t) (Eq [23]), figure [H] 
shows that a tomogram based technique gives a much better agreement with 
the known analytical phase derivative. 

For the three components of the plasma signal we have no way to di- 
rectly verify the accuracy of the tomographic estimates, because the com- 
puted phase derivative is not exempt from noisy corruption. But in any case 
the corresponding spectrogram plots show (Figs{T5l [16] and [17]) that the to- 
mogram allows for a good time resolution and in no case departs from the 
approximate values obtained in the corresponding spectrograms. 

6 Remarks and conclusions 

The tomographic technique for component analysis and computation of the 
phase derivative seems to provide an useful tool for the analysis of refiectom- 
etry signals. The component separation technique contains more information 
than the classical filtering techniques that have been used in the past. In ad- 
dition, the TDM method of phase derivative calculation associated to LM 
filtering compares favorably with those obtained by spectrograms. Neverthe- 
less, a few issues must be addressed: 
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Figure 8: Modulus and real part of the second component of the reflectometry 
signal, corresponding to the reflection on the plasma. For visual puposes, the 
average of the modulus is shifted by 0.2. 

1. How many components should be separated? [9] 

From the tomogram itself one must decide how many components should 
be extracted from the signal. From the analysis of a great number of re- 
flectometry signals it turned out that in some cases the third component, 
corresponding to the multirefiections, was very weak. Maybe, in this case, 
only two components should be extracted. 

2. Separation of the components: for which should the separation be 
performed? 

For 6* ~ 0, the spectrum {cn{9)} is very close to the time representation 
of the signal. Then, the coefficients Cn{0) are almost all different from 0, and 
it is not possible to make the separation of the components. For ^ ~ f, the 
spectrum {c„(6')} is very close to the frequency representation of the signal. 
Then, many coefficients Cn{0) are equal to and the separation can be per- 
formed. But it is not the best choice. The best choice for 6q is where the 
spectrum still has many non null coefficients Cn{0) and where it is still possi- 
ble to make the separation by looking for concentrations of the tomographic 
probability. In the case of the reflectometry signals, the best choices seem to 
be around Oo = % (see Figl5]). 

3. How is the phase derivative dt(j){t) extracted? 

3.1 First, one uses the time representation of the components to decide 
if all parts of the signal are relevant, or if some of it is just flltered noise 
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Figure 9: Phase derivative of the second component of the reflectometry 
signal for t G [11, 20] s. TDM and the TDM+LM filtered phase derivative 
(LMi5 : bold) are presented on the same plot. 

(this is the case for the initial time interval in the second component of the 
reflectometry signals) 

3.2 Then, extract the phase derivative using TDM on the tomogram for 
6*0 (For the first component of the reflectometry signal this was sufficient to 
extract the phase derivative which is almost constant). 

3.3 The use of LM filtering on the phase derivative can be relevant. For 
the second and third component it seems necessary to apply a LMi^ low pass 
filter on the phase derivative. 

4. The refiection on the porthole can be used to calibrate the measure- 
ments. This reflection can be detected after the time ti corresponding to the 
traveling wave from the emitting antenna to the receiver antenna. The group 
delay r^i of the first refiection, computed from the phase derivative should 
be equal to ri. 

The calibration of the measurements can be done by shifting the group 
delay Tg, obtained for each component from the phase derivative, by the 
quantity At = Ti — Tgi. 
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Figure 10: For t G [ll,20]s, LMi^ filtered phase derivative of tlie second 
component of tlie refiectometry signal obtained with TDM (thin line) and 
TNM (bold line). 

given us access to the refiectometry data. 
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Figure 11: Modulus (bold) and real part of the third component of the 
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visual puposes, the average of the modulus is shifted by 0.06. 
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Figure 12: Phase derivative of the third component of the refiectometry 
signal, corresponding to multi reflections on the vessel, using TDM+LM 
(LMi5)on the tomogram for ^ = 
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Figure 13: Phase derivative of the three components of the reflectometry 
signal, obtained with TDM+LM (LM15). The mean value of dt(pi{t) is 
equal to -18 rd/s, dt(p2{t) to 75 rd/s and of dtip^it) to 140 rd/ s. 
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Figure 14: Phase derivative of yR{t) defined by ( !23l) using TDM+TBD on 
the tomogram for 6^ = |. The bold line is the analytical curve of the phase 
derivative. For comparison the dotted lines correspond to the maxima of a 
moving window FFT spectrogram with a resolution ±9.9r(i/s (see section 5). 
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Figure 15: For the first component of the reflectometry signal, corresponding 
to the reflexion on the porthole phase derivative obtained by TDM+LM 
(LM15). The dots are the maximum of a moving window FFT spectrogram 
with a resolution of ±9.9r(i/s . 
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Figure 16: For last part {t G [ll,20]s) of the second component of the 
reflectometry signal, corresponding to the reflexion on the plasma : phase 
derivative obtained by TDM+LM (LM15). The dots are the maximum of 
a moving window FFT spectrogram with a resolution of ±9.9rd/s . 
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Figure 17: For the third component of the reflectometry signal, correspond- 
ing to multi reflexion phase derivative obtained by TDM+LM (LM15). The 
dots are the maximum of a moving window FFT spectrogram with a resolu- 
tion of ±9.9r(i/s. 
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